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We derive the rapidity evolution equation for the gluon four-point function in the dilute regime 
and at small x from the JIMWLK functional equation. We show that beyond leading order in N c 
the mean field (Gaussian) approximation where the four point function is factorized into a product 
of two point functions is violated. We calculate these factorization breaking terms and show that 
they contribute at leading order in N c to correlations of two produced gluons as a function of 
their relative rapidity and azimuthal angle, for generic (rather than back-to-back) angles. Such 
two-particle correlations have been studied experimentally at the BNL-RHIC collider and could be 
scrutinized also for pp (and, in the future, also AA) collisions at the CERN-LHC accelerator. 

I. PRODUCTION OF TWO CORRELATED PARTICLES 



The evolution of QCD amplitudes with energy is described by the Balitsky hierarchy [l[ or, equivalently, by the 
CN| i JIMWLK 0] functional renormalization group equations. They essentially represent generalizations of the well-known 
BFKL equation Q for the evolution of the two-point function to evolution equations for arbitrary n-point functions 
including the non-linear effects due to high gluon density. In the unitarity limit of high parton density the Balitsky 
hierarchy is not closed: the derivative of any n-point function with respect to energy (or rapidity Y ~ log E) involves 
all rra-point functions (m > n). In the dilute regime, however, the hierarchy can be truncated to obtain closed evolution 
$^H, equations for each n-point function. 

Prior work in this field has mostly focused on the evolution of the two-point function and its perturbative unita- 
rization at high energies. The purpose of this paper is to point out that information on the four-point function could 
be obtained from two-particle correlations in inelastic high-energy collisions in a certain kinematic regime (see below) . 
Moreover, we argue that the B- JIMWLK equation for the four point function can not be factorized as a product of 
two BFKL two-point functions. We show that the terms that violate this factorization actually contribute to the 
' correlation function at leading order in N c . 

We consider the correlation of two particles with transverse momenta pj_, q± (we shall drop the subscript _L from 
■ now on to avoid cluttering of notation) and rapidities y p , y q , respectively: 



O C(p q) - / ^ ) ( dN )( dN ) (1) 

\d 2 pdy p d 2 qdy q / \d 2 pdy p / \d 2 qdy q / 

The brackets denote an average over events and the momentum distributions shall be normalized according to 

where (N) is the total average multiplicity per event. It has been argued in ref. Q that in the high-energy limit 
(but fixed p, q, y p , y q ) the leading contribution to C(p, q) is due to diagrams such as the one depicted in fig. [TJ For 
these diagrams the hard amplitudes are disconnected but the correlations arise because for either one (or both) of 
the colliding hadrons the ladders in the amplitude and/or the conjugate amplitude connect to the same color source. 
These two-point functions are essentially the unintegrated gluon distributions of the hadrons; they are of order 1/g 2 
when the transverse momentum in the ladder is below the saturation momentum Q s of the corresponding hadron. 

Diagrams such as fig. Q] should dominate C(p,q) even at high (but not asymptotically high) transverse momentum, 
p.q > Q s , provided one considers generic relative angles cos<p = p ■ q/{\p\ \q\) (in particular, away from the region 
of "back-to-back" jets, <f> ~ 7r). On the other hand, at leading order in a s , when p,q 3> Q s the gluon pair should 
originate from the same ladder; when the rapidity difference between the two produced gluons and the two beams 
are smaller than ~ l/a s the ladder is DGLAP-ordered and C(p,q) should approach 8(p + q) (back-to-back dijet). 
When \y p — y q \ > l/a s the delta- function is smeared out by a BFKL-ordered ladder inbetween the produced gluons 
(Mueller-Navelet jets @). Instead, here we consider the situation where p,q are somewhat larger than but on the 
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FIG. 1: Correlated production of two particles with generic relative azimuthal angle at leading order. The blobs denote 
the unintegrated gluon distribution of the projectile A or target B, respectively, and the light-cone momenta are £1,2 = 
(p/ v / i)exp(±j/ J> ), 2:1,2 = {q/y/s)exp(±y q ). 



order of Q s ; also, \y p — y q \ should be significantly smaller than the total rapidity window between the two beams, 
and the relative azimuthal angle <fi <C 7r. When p and q are on the order of a few GeV it is necessary to subtract the 
background of uncorrelated particle pairs to reveal the structure of the correlation function. 

We note that two-particle correlations away from the back-to-back regime have recently been measured at the 
BNL-RHIC accelerator at ^/s — 200 GeV (per colliding nucleon pair) for proton-proton, deuteron-gold, and gold-gold 
collisions 6-8]. For the former systems only a narrow peak due to fragmentation of the triggered parton have been 
observed. For collisions of heavy ions, on the other hand, C(p,q) exhibits a "ridge" -like structure: it is narrow in <fi 
but extends over several units in Ay = \y p — y q \. The absence of measurable correlations in pp and d + Au collisions 
may be due to the smallness of the saturation momentum Q s for a proton or deuteron at RHIC energy. Also, the 
measurements from RHIC might be expected to be rather sensitive to the initial conditions for the evolution equation 
at moderately small xq. At the higher energies of CERN's LHC collider, the saturation momentum of a proton 
measured from the central rapidity region is expected to be on the order of I GeV and such correlations could be 
sufficiently strong to provide information about the QCD four-point function at small x. 

The diagrams like the one from fig. [T] arise from factorization of the four-point functions in the field of the projec- 
tile/target into products of two-point functions Q (unintegrated gluon distributions). Doing so, however, picks up 
only the leading- jV c contribution to the four-point function. More generally, C(p,q) is given by 

/ dN 2 \ _ g 12 f A d 2 h L ll {pM)L^{pM)L v {qM)L v {qM) 

\<Ppdy p cP q dy q / ~ 64(2tt)« KJavfJtfWlgcc,}***) J (2^ k? (p - (p - krf (q - krf (q - k 4 f 

x (pl(fc2)p*A(fc4W c (fciW d (fc 3 )p*B(p- k 2 )p* b B (q- k 4 )p B c '(p- k x )p B d '{q-k^)) (4) 

~~ 64(2^)6 Khaa'Jg'Wjgcc'Jg'dd') J JQ ^ _ _ ^ _ _ ^2 

x (p'^Wa^W^iW^)) (p*b(P~ k2)p*%- k 4 ) PB c '(p~ h)p B d '(q- hajft) 

In the second step we have assumed factorization of the wave functions of projectile and target. denotes the 
Lipatov vertex which satisfies 

L„(p, k^L^p, fe) = [S t3 S nm + e«e nm ] k\(p - k,y k%(p - k 2 ) m (6) 
4fc 2 

L^p,k)L»(p,k) = r ( P -fc) 2 • (7) 

The expression ([5]) is depicted in fig. [5J Here, p(r) denotes the color charge density per unit transverse area at a 
transverse coordinate r and p(k) is its Fourier transform. Its two-point function is related to the unintegrated gluon 
distribution &(x, k 2 ) via 

(p* a {k)p\k')) (x) = ±^—(2nfS(k k') $(x, fc 2 ) . (8) 



FIG. 2: Correlated production of two particles with generic relative azimuthal angle at leading order. The blobs denote the 
four-point functions for the projectile A or target B, respectively. 



With this normalization one recovers the LO fc^-factorization formula for the single-inclusive distribution from the 
diagram [3] with the standard prefactor [9(: 



dN 
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(9) 



where cto is the transverse area of the collision (note that in our convention §(x, k 2 ) is the density of gluons per unit 
transverse area and it therefore contains a factor of 1/<tq). 
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FIG. 3: Single-particle production from fcx-factorization at leading order. The blobs denote the unintegrated gluon distribution 
of the projectile A or target B, respectively. 

In a mean held (and large N c ) approximation one may factorize the four-point functions from eq. §5§ into products 
of two-point functions, 

(p a p b p c p d ) - S ab S cd {p 2 ) 2 + S ac S bd (p 2 ) 2 + 5 ad 5 bc {p 2 ) 2 + • • • , (10) 

where p 2 = (pp), and the momentum dependence of the two-point function has been suppressed. Then, one of the 
nine contractions corresponds to the square of the single-inclusive distribution: contract the first p with the third and 
the second with the fourth, for both projectile and target. The color factor for this diagram is 1 

fgaa'fg'bb'fgcc'fg'dd' (p*APA C ) (p^aPA^ (p*°bPB C ^ (p^PB^ ) (11) 
~ fgaa>fg'wfgcc'fg'dd'5 aC 5 bd 5 a ' C '5 b ' d ' = N 2 (N 2 - if . (12) 



1 Not including factors of N c which will enter once (pp) is expressed through $ via eq. JHJ. 
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The remaining eight diagrams correspond to a color factor of (we take fig. Q] as an example) 

fgaa'Uwfgco'Udd' (p* A P A C ) (p^aPA 4 ) ( P *% P*b) (pb" W) (13) 

~ f a aa>U bb >f g cc>U dd '6 ac S bd S a '»'5 c ' d ' = N 2 (N 2 - 1) . (14) 

Thus, two-particle correlations are suppressed by a factor of N 2 — 1 as compared to uncorrelated production. For this 
reason, the leading- iV c ansatz (|10p may not capture the complete result for C (p, q) . Below, we derive the evolution 
equation for the four-point function from JIMWLK. We determine the corrections beyond the mean-field and large- iV c 
approximations to the rhs of ([TU|) and show that these corrections contribute at the same order in N c to the correlation 
function. 

In this regard, wc should point out that N c corrections to the two-point function in the dense regime were found to 
be exceptionally small [10]. However, this needs not be true for the four-point function. In fact, we shall argue below 
that we do not expect N c corrections to the four-point function to be anomalously small, even in the dilute regime. 
A verification or falsification of this expectation via exact numerical solutions would be very valuable. 



II. EVOLUTION EQUATION FOR THE FOUR-POINT FUNCTION 

In this section we present the equation describing the rapidity evolution of the four-point function (a^apa^ag) 
obtained from the JIMWLK equations, which include terms of subleading order in N c . In this context it is more 
natural to work in coordinate space, so r, s, f, s denote transverse coordinates; the four-point function in momentum 
space can be obtained by Fourier transform. We also find it preferable to work with the fields a rather than the color 
charge densities p; at leading order and in covariant gauge, they are related in coordinate space by 



A^{x + ,r) = S tl -a(x + ,r) = -g 5^ S(x+)— T p{x + ,r) , 



vi 



(15) 



for a hadron moving at the speed of light in the negative z-direction. Since this field also satisfies A + — 0, the only 
non-vanishing field-strength is F~ l — —d l a. In momentum space we have the relation k 2 a(k) = gp(k). 

The JIMWLK evolution equation for the four-point function to lowest order in the fields can be shown to be (see 
appendix) 



dY 



/a b _ c _ <i\ 



7T 



d 2 Z 



(2tt) 3 



d 2 : 



(r — z) 2 
z) • (f — z) 



a a r a h z a c s ai 
(f — z) 2 



(s-z) 2 



(r — z) 2 (r — z) 2 



I rena rfnc ( r Z ) ' (f Z ) r e / 

1 J (r - z) 2 {s - z) 2 [ r " 



+r a f 



(r 



(5 



e f 

a r a J s 



(r — z) 2 (s — z) 2 

z ) ■ (S - z) r e / 



(r — z) 2 (s — z) 2 



, feKbjfKd (r-z)-(s-z) 
(f — z) 2 (s — z) 2 



+ r c f 



fKd ( S ~ Z ) ■ ( S ~ Z ) 

(s — z) 2 (s — z) 2 



J 



U 2 U S 



a e z a f s 



a e z a{ 



a,?a,% 



{s-z) 2 



a z a{ 



<4 



ot e z a{ 



a b F a c s 



a e z a f z ] a a r a d s 



a.a; 



0&0& 



a a r a b F a c s af 



(16) 



This expression neglects contributions from higher n-point functions on the rhs; in the dilute regime, i.e. when the 
transverse momenta of the produced particles are higher than the saturation momenta of the colliding hadrons, this 
approximation should be justified. 

In order to derive the color structure of corrections beyond the large- N c approximation, we factorize the product of 
four point functions on the rhs of eq. (|16l) into products of two point functions. This Gaussian approximation reduces 
the evolution equation for the four point function to a product of two BFKL equations (for the two point function) 
plus extra terms which provide corrections to the factorization (|T(J|) . The result is 



dY 
where 



(a a r a b a c s a d ) 



_d_ 

dY 



[S ac 5 b 



5 ab S cd a 2 ,a 2 



5 ad S bc a 2 



-,a 2 _ 



1 as I j2 f T^abcd I j^abcd . 



(17) 



jpabcd — 



abed j-anb rctzd 

jarzb jeftd 
jand jbnc 
janb jend 



(r - s) 2 



jand jb 



(r 



(r - s) 2 



(r — z) 2 (r — z) 2 (s — z) 2 
1 (f - r) 2 



2(r 

2 



1 



-z) 2 (f-z) 2 2(r 
g + (r O s, s 4-> f) 
(s-f) 2 



+ (r o s, s o r) 



rpabed 



(r — z) 2 
(r- 



(r 



z) 2 (r — z) 2 
1 



(s-z) 2 



(r — z) 2 (s — z) 2 (r — z) 2 (s — zf 



(r — z) 2 (s — z) 2 



a 2 a 2 



+ (r O s, s O r) 
(r O f) 



In (|T8|) all terms in Fq, F\ and F2 are to be duplicated with the substitutions indicated explicitly in the brackets. 
Second, all terms in Fq and Fi are to be duplicated again substituting a <h> b and r -H- f , Then, all terms in Fq and 
F 1 should be duplicated a third time exchanging c O d and s ■(-> s. Furthermore, all terms in F 2 are to be duplicated 
while letting b c and r O s. Finally, the terms obtained in the last substitution (only) in Fi should be duplicated 
exchanging c O d and r ■<-> f. 

The first term in (jTTJ) provides the leading-^ contribution to the four-point function. The second term gives 
corrections beyond the large- N c factorization (|10p . Since an analytic solution to the evolution equation for the four 
point function is not within our reach, a numerical investigation of these terms and their magnitude would be extremely 
useful. Nevertheless, from 



d Y (p a p b p c p d ) - a s N c S ab S cd (p 2 ) 2 + a s f acK f bdK (p 2 ) 2 , with p 2 (Y) - e Q ^ Y 
one might expect that, generically, the solution to this equation has the following color structure: 
(p a p b p c p d ) - S ab 6 cd {p 2 ) 2 + 8 ac 8 bd {p 2 ) 2 + S ad S bc {p 2 ) 2 + 



(19) 

(20) 
(21) 



(Note that the various two-point functions depend on different coordinates/momenta and so each of the above terms 
is distinct.) The color factors emerging from the products of the Kronecker tensors have already been discussed above, 
eqs. (fT2"j) and (|14[) . However, some of the products of a leading- N c term from the first line ([20)) with a subleading-7V c 
term from the second line (fl?Tj) also contribute at the same order N 2 (N 2 — 1). For example, 



1 



ra c rb d' _c r e r rabn rcdn 

J gaa' J g'bb' J gec' J g' dd' J j 



abK fcdK 



N c 5 ac S oa f a0K f 



N 2 (N 2 -1) 



(22) 



The other products are worked out in appendix [C] 

This shows that some of the subleading-iV c contributions from the four-point function actually enter C(p,q) at 
leading order, compare to eq. (TH| . Previous results from the literature [H (also see [TTJ| ) are therefore not complete. 
Nevertheless, the correlations described here should still extend over several units in \y p — y q \ [ll|. Quantitative results 
for the JIMWLK four-point function and for the corresponding two-particle correlations C(p, q) as functions of the 
transverse momenta p, q, relative azimuth <f> and relative rapidity \y p — y q \ remain to be found. 

In summary, we have argued that two-particle correlations from high-energy collisions may provide some insight into 
the QCD four-point function. This should be the case, in particular, when the transverse momenta of the produced 
particles are not very much higher than the saturation momenta of the colliding hadrons and when their relative 
azimuthal angle is sufficiently less than tt. The narrow (in both azimuthal and polar angle) jet-like fragmentation 
peak should sit on top of a "background" which is broader in the relative rapidity \y p — y q \. 

If expanded in powers of N c , the leading contribution to the four-point function is given by the product of two 
BFKL two-point functions. However, we find that genuine B- JIMWLK subleading-iV c corrections also appear in the 
correlation function C(p, q), at leading non- vanishing order in N c . The correlations mentioned here represent an 
interesting opportunity to study the non-trivial structure of the four-point function of the B- JIMWLK hierarchy, 
both theoretically and experimentally. 
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Appendix A: Two-point function and BFKL 



The JIMWLK equation for the two-point function is 

4?KO = \l d 2 xd 2 y4jrrt d v 4^a> 



dY 



Sa b lxv Sa d r s 



—n bd 
5ai lxy 



where 



bd 

lxy 



d 2 xd 2 y 
1 f d 2 z [x — z) ■ (y — z) 



5 a c 1 

a a + - 

Sa d r s 2 



d 2 xd 2 vn bd — —a a a c 



7T J (2tt) 2 {x~z) 2 {y-z) 2 
We start with the first term from eq. (IA2[) : 



Sa b 



bd 

lxy 



l 



d 2 z [x — z) ■ (y — z) 



ir J (2tt) 2 (x-z) 2 {y-z) 2 
-ig f d 2 z (x — z) ■ (y — z) 



V *lk V *- V *lk V *- V 'lk V « 



7T J (2tt) 2 (x-z) 2 {y-z) 2 
Note that in (IA4I) the derivative does not act on the on the left because 
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[S(x - y)V^V y t b - 8(x - z)V}V z t b - 5{x - y)V}V v t l 



, bd 



d^M x ) bd = ig5[x ~ u) {tbylx)bd ~ {tb)be {vlx) * d = ° 



(Al) 
(A2) 

(A3) 

(A4) 
(A5) 

(A6) 



The first term in the bracket from eq. (| A5|) vanishes because S(x — y)V^V y t b = S(x — y)t b and (t b )bd = 0. To apply 
the same argument to the second term from (|A5|) we rewrite 

2(x - z ).(y- z ) = (x- z) 2 + (y - z) 2 - [(x - z) - (y - z)] 2 = (x - z) 2 + (y - z) 2 - (x - y) 2 . (A7) 

Hence, the second term from (|A5j) becomes 



d 2 z (x — z) 2 + (y — z) 2 — (x — yf 



8{x-z) [VjV z t b ] M = . 



(A8) 



2tt J (2tt) 2 {x-z) 2 (y~z) 2 

The first term again vanishes when x = z while for the other two terms the divergences at x = z cancel and so they 
vanish due to the color structure. 

We can therefore simplify eq. (|A5I) to 

d 2 z (x — z) ■ (y — z) 



For the second term from eq. (|A2I) we need 



(2tt) 2 (x-z) 2 (y-zf 



[V}V y t l 



, bd 



= 8 ad (r-y)a<; + 5 cd (s-y)a a r 
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x y 



(A9) 

(A10) 
(All) 



Eq. (|A2I) turns into 
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To expand the rhs to second order in the fields we need the following expressions: 



i ba 
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(A17) 
(A18) 
(A19) 
(A20) 
(A21) 

(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



In the step from (|A23j) to (|A24j) we have used the rapidity ordering of the fields; hence, only one of the ~ a 2 terms 
can contribute. 

Using these expressions in (|A16[) gives 



dY K r s * (2tt) 3 J {r-z) 2 
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We now take the expectation value on the rhs using 



/ a c\ eac 2 



(A27) 



(A28) 



where a 2 _ s is essentially the unintegrated gluon distribution function. This turns eq. (1A27|) into (we drop an overall 

8 ac ) 



dY Ur - s 



2g 2 N c 
' (2tt) 3 



d 2 z 



r- 



-\a 2 ~a 2 1 

2 I T— 8 ^Z — Sj 



2g 2 N c 
(2tt) 3 



(s — z) ■ (r — z) 
(s — z) 2 (r — z) 2 



To cast this into the familiar form for the BFKL equation we symmetrize the kernel of the first term and rewrite that 
of the second term in a different form: 
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Then, the part of (|A29|) involving a 2 _ s (the BFKL virtual part) becomes 

g 2 N c f , 2 (r-s) 2 2 a f , 2 (r-s) 2 



(2tt) 2 
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(A30) 
(A31) 

(A32) 



where a = a s N c /n. 

The other terms from (|A29[) where a 2 depends on the integration variable z (the BFKL real emissions) turn into 
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Appendix B: Four-point function 



The JIMWLK equation for the four-point function is 

-^p{a a r a b a c s a d ) = ~J d 2 xd 2 y^r] e J y jja a r .a b f a c s ai 
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5 2 
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(Bl) 
(B2) 



We begin with the first term. The square bracket can be taken over from eq. (|A9|) . and may also be written in the 
form 
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The first term on the rhs of (|B2|) becomes 
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To linear order in the fields, 
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and so eq. (|B5I) is equal to 
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(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



This is the final result for the first term from (|B2|) . In the next section [Bj] we will simplify this further by taking the 
expectation value with a Gaussian weight. 

We now turn to the second term from eq. (|B2p . With rj^ b — r) b ® we find 



d 2 xd 2 y rfl T a a r a b a%a d = tf b a c s a d + rj%a$<4 + rf r d a\a c s + 4»t + 4t a > C s + ■ ( B9 ) 



xy Sa% Say 



The expansion of r/^ b in powers of the field starts at second order. From (|A3|) . 
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We note that the two terms a 2 exhibit different ordering in rapidity as they arise from the expansion of V z and V% 
to order g 2 , respectively. For what follows, we combine them into a single ~ a 2 term which shows up no matter how 
rapidities are ordered. Using (jBllj) in (|B9|) we obtain 
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This is the final result for the second term from (|B2[) . In the next section we compute its expectation value for a 
Gaussian weight. 



1. Gaussian approximation 



To simplify the expressions further, we assume that the expectation value of the four-point function on the rhs of 
the evolution equation is taken with a Gaussian weight so that it is given by a sum over all possible Wick contractions: 
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where we used (|A28j) . The above factorization into two-point functions reproduces the evolution of the four-point 
function at leading order in N c . The lhs of (IB2[) then becomes 
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We continue to analyze only one color channel corresponding to the first term ~ S ac 5 bd (the others are similar). The 
two-point functions satisfy the BFKL equation and therefore 
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Thus, (|B15|) turns into 
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At leading order in N c the evolution of the four-point function should be determined by the BFKL evolution of the 
two-point function, as given in the previous equation [plus similar terms from eqs. (|B16IB17[) ]. In what follows, we 
verify this explicitly from the full JIMWLK evolution equation, eq. (|B8|) plus (|B12[) . However, we also derive the 
contributions which are subleading in N c . 



a. Virtual terms from JIMWLK 



Forming pairwise contractions on eq. (|B8[) gives 
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All of these terms will cancel against corresponding pieces from eq. (|B12[) . 

We organize the virtual contributions (where none of the fields depends on the integration variable z) from eq. 
according to the various possible contractions. The term involving ol 2 _ s ol!^__s is 
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The first term on the first line cancels the corresponding term involving a 2 _ s a 2 _ s from eq. (|B21[) . The rest of that 
line contributes to the rapidity evolution of the four-point function at order 0(N C ) and is seen to match the first 
line of eq. (|B20|) , i.e. the "virtual correction" part of the leading- N c BFKL equation for the four-point function. The 
second line in (|B24[) is a correction which is suppressed by a relative factor of 1/N C . 
Along the same lines one finds that the contribution from eq. (|B12|) involving 
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The first term on the first line will again cancel with a similar term <- 
eq. (HU. 

Finally, we also list the contribution from eq. (|B12[) involving a 2 _ g a 2 _ 
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As before, the very first term in this expression will cancel. To summarize the results obtained so far: the evolution 
equation for the four-point function in mean-field approximation is given by 



d 

dY 



I „ a b c d\ 

[a r a f ,a s a s ) 



dY 

9 2 



(S ac S l 



i a 2 -a 2 



, cab cca A A , caa cbc A A \ 

f-s-r-s + od a s _ s a r _ f + 6 6 a f _ s a r _^) 



(B27) 



(2tt) 3 

jclkc jbnd 



d 2 z< 



(r — r) 



(r-s) 



(s-ff 



(s-sf 



_(r 


— z) 2 (r — z) 2 


(r 


-z) 2 (s-z) 2 


(s 


— z) 2 (r — 




(s 


-z) 2 (s- 


z) 2 _ 




(r-s) 2 




(r-s) 2 




(s-f) 2 






(s-f) 2 




_(r 


— z) 2 (s — z) 2 


(r 


-z) 2 (s-z) 2 


(s 


— z) 2 (f — 


z) 2 + 


(5 


— z) 2 (f — 


z) 2 _ 




(r — f) 2 




(r - s) 2 




(s-f) 2 






(s-sf 




[(r 


— z) 2 (f — z) 2 


(r 


— z) 2 (s — z) 2 


(3 


-z) 2 (f- 


z) z 


(s 


-z) 2 (s- 


z) 2 _ 



2 2 

x a- - 

r—sr—s 



2 2 
^■r — r s—s 



2 2 
r — s^r— s 



The first line, eq. (|B27[) , involves only standard BFKL evolution of the two-point functions and provides the leading- 
N c contributions to both real emissions and virtual corrections. These terms can be written in explicit form by using 
the BFKL equation for the rapidity evolution for the two-point functions; for example, the first term from that line 
is given in (|B20I) . The remaining lines correspond to the subleading (in N c ) contributions to the virtual terms for 
JIMWLK evolution. In the next section, we derive the corresponding real emission terms. 



b. Real terms from JIMWLK 

We begin with terms which contain only one z-dependent field, as given by the second and third columns in (|B12[) . 
The first line is given by 
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The second plus third terms from the first line of (|B12[) thus becomes 
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The remaining terms from ()B12|) with one z-dependent field are easily obtained via permutations of color indices and 
coordinates. The full result is 
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Adding now the terms from eq. (|B21[) with one z-dependent field, and those from eq. (|B12|) with two z-depcndcnt 
fields, leads to the sum of the following two contributions; the terms proportional to N c are given by 
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These provide the leading-iV c contributions and correspond to BFKL evolution of the product of two two-point 
functions; compare to eq. (1B20[) . The genuine JIMWLK contribution is 
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This expression (|B34[) . which is to be added to the rhs of (|B27|) . is our final result for the evolution equation of the 
JIMWLK four-point function for a Gaussian weight. It can be rewritten in more compact form by exploiting some of 
the symmetries of (a^a b a^a d ) , eqs. (|17|18|l . 



Appendix C: Color factors 



In this appendix, we compute the color factors for the products of a leading- iV c term from the first line (f2"U|) with 
a subleading-A'c term from the second line (|21|) ; color indices from the target side carry a prime and we need to also 
include the remaining structure constants from eq. ([5]). Using the following SU(iV c ) identities 

(CI) 
(C2) 

(C3) 

(C4) 

(C5) 
(C6) 

■^6 a '»'6 c ' d ' {f gaal f gcd r K ) (f g 'wf g ' dd >f bdK ) = \n 2 (N 2 - 1) , (C7) 





jabK jcbn 


= N c S ac , 




jabn jccLk 






d acK d acr 


= ( n '-k) 5 "- 


jgaa jgcc jacK, 


= i 3 tr adj < a W 


c £CL C K 

2 J 


ra rc fd fb 
J a' gJ gc' J c' g' J g' 


a> = tr a dj t a t c t d t h 


= S ac S bd + S ab 5 cd + — (d acr d dbr - d adr d bcr + d abr d cdr 




d adK d bcK d acr d bdr 


= 0{Nt) , 



one derives 



ra b rc d £ £ £ £ £abhz £ 

J gaa' J gcc' J g'bb' J g' dd' J J 



N, 



abhz £cdK 



N c - 



f afcK / cdK tr ad ,t a t c t d t b 



jabn jcdn | ^ac^bd _|_ ^ab ^cd _|_ ^ c ^acr ^dbr ^adr ^bcr _|_ ^abr ^cdr s 



gbd _|_ _ jabtz jcdK ^acr ^jbdr 



N 



2 V N r 



0{Nl) 



(N 2 C - 1) (C8) 



— S a ' b 'S C ' d ' 'f g aa>f g cc'f g 'wf g ' dd 'f adK f bCK = J^f adK f bCR tr adj t a t C t d t» 

N 



= _fad K jbc K ( S ac S b d + s ab S c d + _£ (jacr jdbr _ jadr jbcr + jabr jcdrj j = q 



1 

N~ c 
1 

N 
l 

nZ 



5 a ' c '6 b ' d 'f gaa ,f g , w f gcc ,f g , dd , /««/««« = N c S ac S oa f acK f 
5 a ' c '5 b ' d ' f 9 aa>f g iwf 9 cc>f g > dd > f abK f cdK = Nc s^s bd f abK f dK = N 2 (N 2 - 1) 



zac sbd fcicK £bd-K 



ra c xb' d £ £ £ £ £ddK £Ock at cac cod £.aaK £ 

J gaa' J g'bb' Jgcc' J g' dd' J J = ^ C J J 



adK £bcK 
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N 2 (N 2 - 1) 



N, 



1 8 a ' d '8 b ' c \f gaa ,f gcc> r*)(f g , bb ,f g , dd ,f bd *) = -jN 2 (N 2 -l), 



_ SCl'd' $tb' c' £ £ £ £ i-iihi: r 

j\j J gaa' J g'bb' J gcc' J g' dd' j J 



abn. £cdn 



N, 



ca = o 



^5 a ' d 'S b ' c 'f gaa ,f g , bb ,f gcc ,f g , dd ,r d *f b ™ = -±N 2 (N 2 -1). 



(C9) 
(CIO) 
(Cll) 
(C12) 
(C13) 

(C14) 
(C15) 
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